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Abstract 

The issue of how superconductivity originate in the Cu02 planes believed to be crucial to un- 
derstanding the high Tc superconducting cuprates is still an going debate. In the wake of recent 
experimental observations of the the Zhang-Rice singlet (ZRS), its formation and propagation need 
to be revisited especially by using a simple approach almost at a phenomenological level. Within 
a highly simplified correlated variational approach (HSCVA) in this paper, a new formation of the 
ZRS as constituting the ground state of a single-band t-J model of the Cu02 planes is developed. 
This formation is then used to demonstrate how the ZRS can be propagated as a probable Cooper 
channel in the Cu02 planes. 



PACS numbers: 74.72.-h, 74.25.Jb, 74.20.-z 



I. INTRODUCTION 



One early consensus after the discovery of the high Tc superconducting cuprates, is that 
the key to understanding these materials is the Cu02 planes common to all of them. How- 
ever, after two decades of intense study, there is no consensus understanding yet of the origin 
of superconductivity from this feature. Interestingly, it was experimentally demonstrated 
early enough that when the parent material is doped, the mobile holes reside on the O site 
of this plane. This led to the suggestion that a three-band Hubbard model (i^s^) in which 
;he hole is mobile and carry a spin should be the starting point to investigate these material 
l|. In a seminal paper j2|, Zhang and Rice suggested that this hole on the O site will form 
a singlet state with a hole on the Cu site at the central of each Cu02 plaquette to form a 
single-band character. This is the mapping of the H^b into the single-band Hubbard model 
(/fife) when J = At'^/U which is consistent with the Anderson's proposal that a strong on-site 
Coulomb interaction among a partially filled band of Cu 3d levels should be the starting 



effective single-band model 
follow this line of thinking 



:or_the supeconducting cuprates. Since then, the researchers who 
3, 



, J] believe the ground states of the Cu02 are the Zhang-Rice 
singlet (ZRS) which are expected to become the Cooper pairs of the superconducting states 
when liberated from the insulating host material. However, there have been opposing views 
on the equivalence of the ifat and Hn, and also that the added holes do not hybridize into a 
ZRS jsL |6|. These views may need some reconsideration in the wake of recent experiments 



n 



[7, ISl, l9|] and calculations [10|, lUl, ll2| which have demonstrated the existence of the ZRS. 
Meanwhile among the workers starting from the ZRS, there is still no consensus on the 
appropriate parameters to add to the t-J model to obtain an effective Hamiltonian for the 



superconducting cuprates 



13 



14| . The reason is that the effects of these parameters even in 



first principle calculations, depend on the starting CuQ2 configuration. It follows then that 
a simplified determination of the possible configuration at a phenomenological level but not 
bias to experimental results can give useful insight into the starting structure of the Cu02- 
This is the purpose of this paper. 



Interestingly the ZRS is obviously a two-hole state [4, Il5|, llGj and this makes it straightfor- 
ward for us to apply to the CUO2 planes our recently developed formulation of the statistical 
equivalent of the Hubbard (t-U and t — t' — U) model using a highly simplified correlated 
variational approach (HSCVA)[17]. The simplicity of this approach makes it easy to give 
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useful insight on how the competition between itineracy and locahzation in strongly corre- 
lated systems could lead to the exciting properties in these systems in d = 1,2,3. In the 
present formulation, it will be shown how the competition between the motion of a single- 
hole in the the Cu02 planes and its hybridization with the Cu can lead to a Hn^ model. In 
particular, the approach emphasizes the ZRS as the ground states of the CUO2 planes and 
that it can be propagated as a probable Cooper channel in the CUO2 structure within the 
t-J Hamiltonian for the supeconducting cuprates. 



II. METHOD AND FORMULATION 



It is generally believed that the undoped superconducting cuprates is a Mott insulator 
at half-filling (one electron per site). In this state, the system is more stable because the 
electron are reluctant to hop in order to avoid the price of occupying the same site (i.e U). 
Thus the starting CUO2 is viewed as a valence state of Cu "id^ and O 2P^ with a hole at 
the Cu site. One way to induce electronic motion is to create an electronic vacancy on the 
O site which is often called a hole. This free hole attracts an electron of another O which in 
turn creates a vacancy in its place and this new vacancy will attract another electron and 
so on. Thus by this first electron removal (FER) approach, the electron can hop from one 
O site to another, but in doing so, the hole is also hopping from one point to another and 
this is often referred to as the hole motion. As stated above, the formulation in |2l] assumes 
that the ZRS is formed from the hybridization of a hole on each square of O atoms to the 
central Cu to form the local singlet. It follows that this FER approach to investigate the 
propagation of a single hole in the CUO2 plane is a natural starting point to investigate the 
superconducting cuprates and the ZRS states are believed to be the theoretical construction 
of the FER states in photoemission experiments [l5|. 

In the theoretical construction here, we start from the H^b scenario by assuming all four 
possible Cu-0 hybridized states in a plaquette. Since only one hole is present in the FER 
approach, it is argued that only one of the hybridized states will be present at any given time 
to form an effective singlet of each plaquette. Therefore this effective singlet retain the ZRS 
nature of the Cu02 structure in the FER approach. The non-triviality of this consideration 



to the formation here will 



DC made clear below. 



Following the steps in 17|, the basis states of the Cu02 plaquatte in a three-band config- 



3 



Y- 



(a) (b) 




(c) 

FIG. 1: (a) The 2D lattice of Cu-0 layer of hight Tc cuprtates showing the Cu02 planes (b) The 
Cu02 plaquette as a ID ring of O around the Cu so that there is no edge effect hence have periodic 
boundaries (c) Hole doping leading to propagation of the ZRS in the Cu02 planes. 

uration from Fig. lb are , \CuxCr,Oxid' >, \Cu3;Cr,0x20' >, \Cuya,Oyj^a >, \Cuya,Oy.^a >, 
where the a and a represent the spins defined by (t{(t) =1 (|),J, (t)- In numerical index, 
a state can be represented in general with the first site for Cu and the second site for O 
as \iia,jka > where (ii) means there is a particle at the site i that can possibly be in the 
(x,y) direction. Then the state \CuxCr,Oxj^a > becomes \iOa,iiOa > where i = 0, 1, 2, 3, 
...N is the index number of Cu sites. It is worthy to note that the notation used here in the 
configuration of the hybridized states is for convenience of the approach and do not differ 
from the common ones in the literature. Further more, this configuration combined with the 
periodicity of the Cu02 structure implies the eigenstates of the Cu02 plane are eigenstates 
of the translational operator, Tr\'4^ >— exp{lk.r)\ip >, and this gives a recipe to construct 
translationally invariant basis states with additional fixed quantum number, K — ^ which 
will become important to extend the simple approach here in future communication. 
The Cu02 plane is a quasi two-dimensional (2D) structure and therefore requires the use 
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of the 2D kinetic operator, ht, which for a general state, \iia^ka > say, yields the excited 



states simply by adding and subtracting one appropriately 



ht\ii<J,jka >= —t[\{i ± l)ia,jka > ± l)(T,jka > + 

\iia, {j ± l)ka > +\iia,j{k ± l)a >], (1) 

where -t is the nearest neigbour (NN) hopping matrix. 



It has been shown 
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that the variational ground state energy 



for the model Hamiltonian B. = ht (i.e Eq. 1) will yield a matrix 

[ESL'.L + iT)L.L'] = (3) 

where E = Eg/t is the total ground state energy spectrum which can be obtained by 
diagonalization of Eq. (3) and its smallest value is the ground state energy of the systems. 
The Lc in Eq. (3) is the separation between the sites in a basis state written in compact 
form, Lc = \i L' is the new separation of the excited state obtained and T is the number 
of such excited states with V for an operation on any of the basis states. We emphasize 



a, 



the separation depends only on 



that since the Cu-0 distances in the plane are the same 
the direction. 

It is easy to show by diagonalising Eq. (3) for a Hubbard square lattice that the ground 
state energy is the same as the bandwidth W = -8 eV = 2zt (for t = -1) where z = 4 is 
the coordination number. However, Eq. (1) needs to be modified for the hole motion in 
the Cu02 plane. For as it is obvious in Fig lb, the motion is periodic in ID since the O 
sites form a ring around the Cu. An immediate question is why not use the ID hopperer? 
This is not possible because the basis states are 2D. Thus we see right from the beginning 



of this formation why the quasi 2D nature of the CuO_2. 
superconducting cuprates to be smaller than W = 2zt 



2_P 
10, 



ane make the bandwidth of the 



111- 



Eq. (1) for the ID hole motion in the Cu02 plane can be expressed as 

htpp\ncr,jka >= -tpp[\iOa,i(^o±i)Oa > +\Oia,Oi(o±i)^ >], (4) 
where Lc = and L' = ± 1 and tpp is for 0-0 hopping. 

Then taking Eq (4) into account in Eq. (2) will result to a 4 x 4 matrix in Eq. (3) which when 



diagonalisd will yield W = -2 eV. This means the bandwidth of the non-interacting case of 
the Cu02 plane is 75% smaller than that of a normal Hubbard square. This narrow band 
is a key feature in the high superconducting cuprates because it is crucial in determining 
all the possible hopping and interacting parameters that will contribute to the properties of 



10|,ll8|. 



these materials and hence the model Hamiltonian 
Now if we take into account the assumption made above to retain the ZRS in the FER 
approach, Eq. (4) can be expressed in second quantization language as 



pd 



-t 



pd 



Yl Yl K'-Pt-^^r^Pl^ra + H.C 

'- {i} <j,k>e{i} 



(5) 



where d^{d) is the creation (annihilation) of a hole at the Cu Sdx"^ — y'^ orbital and p^{p) is 
the creation (annihilation) of a hole at the O 2p^. and O 2py orbitals, H.C. is the Hermitian 
conjugate while {i} denotes [i,r,r' ,a,a] with r =(x,y) and r' = {y,x). 
It is seen immediately that Ht^^ is the kinetic operator of the t-J model [2I, 
The beauty of the HSCVA is that one can start the investigation by obtaining the band- 
width of the non-interacting case. Thereafter one can introduce the desired interactions and 
monitor their effects. As stated above, the on-site Coulomb interaction of of the Cu holes is 
considered in ^\ to be very large, Udd 00. They therefore considered only the exchange 
interaction between the Cu and O with exchange matrix, J^p. The physical implication of 
this interaction is that as the singlet state propagates, there is also an exchange of spins 
between the Cu and O. Thus a spin exchange distributor, hg which add and subtract one 
spin appropriately to the spins of the basis states will lead to new states which are just the 
antisymmetric states of the basis states. 



hs\iia,jka >= —Jdp\ii<J,jka > 



where 



h., = J, 



dp 



'- {i} <j,k>t{i} 



(6) 



(7) 



Eq. (7) is the XY limit of the anisotropic Heisenberg exchange interaction in second quanti- 
zation language 20!], and it is the interacting part of the t-J model in the CUO2 plane with 
ZRS. Taking Eq. (7) into account in Eq. (2), Eq. (3) becomes 



[EsLcL' + TlcL' + JLcL'] = 0. 



(8) 



where J = Jpdgs/tpdQt with gt and Qs being the Gutzwiller renormahzation factors for the 
kinetic operator and superexchange interactions. The inclusion of these factors is to account 
'or the effect of strong correlations. For in deriving Eqs. (5) and (7), the arguement in 
2l| has to be invoked that since the Cu on-site Coulombic interaction is strong, Udd — ^ 



oo, we will have projected fermions di^ — > Xi^a = (1 — ^4^5.) which do not obey the 
usual commutation rule for free fermions. However, taking into account the Gutzwiller 
approximation, the renormarlized fermions becomes Xida = gdia where g = gt = 2x / {1 + x) 
and g = gs = 4:/ with x denoting the level of concentration of holes. Thus the effective 
Hamiltonian is a combination of Eqs. (5) and (7) with the inclusion of the renormahzation 
factors, 

H^ff = gtHtpd + gshs (9) 

III. DISCUSSION AND CONCLUSION 

To observe the role of the exchange interactions in comparison with Eq. (3), Eq. (8) is 
diagonalized at various values of tpd and Jdp without the strong correlations effect for now; 



it is observed that the exchange interactions has a band narrowing effect 



ll| . For example. 



the common values in the literature [6!], tp^ = 0.4 eV and Jdp = 0.12 eV, gives W = -1.70 eV 



22 1 can be obtained from 



while the experimental value of about 1 eV for the CUO2 structure 
the unrealistic value J dp /tpd = 1.0. However, the inclusion of the strong correlations effect 
by appropriate doping can be used to correct this defect as ~ 1 eV for x = 0.4 for the 
common values {Jdp/ tpd = 0.3) and W ^ 1 eY ioi x = 0.3 for Jdp/ tpd = 0.2 . This is quite 
encouraging as it has been earlier demostrated that the Gutzwiller projected d-wave pairing 



state is most stable for a; < 0.4 |t21j] w^ 
only stable up to x = 0.3 generically 



lich is a fair generalization since that of the cuprates is 
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24l |. The only exception is the Y-Ba-Cu-0 samples 
which has an additional ID Cu-0 chains that plays a role but is not taken into account in 
most calculations in 2D model. Thus the emphasis here is that our formulation has shown 
that appropraite doping of correctly chosen CUO2 configuration within the t-J model may 
lead to results in fair agreement with experiments. It is important to point out that over the 
years, several sophisticated techniques have been used to investigate the inclusion of addi- 
tional hopping and interaction parameters beyond the t-J model to achieve the experimental 
value. Though these studies have contributed to enhanced understanding, there is still no 
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consensus on what parameters should constitute the effective Hamiltonian to study the high 
Tc superconducting cuprates. As stated above, the reason may be due to the starting Cu02 
configuration. The most probable Cu02 configuration is therefore analyzed here. 
It is assumed that the hole is created at the O site Y_ as shown in Fig. (Ic) where the 
thick spin are before exchange interactions and the dotted spin are after the exchange inter- 
actions. Therefore two hybridized dotted spins indicates a ZRS and hence is used here as 
indicative of the ZRS propagation. From the kinetic part of Eq. (5), the hole from the 
will hop to X_ position of the O site between Plaquettes 1 and 2. It can then form a ZRS 
with the Cu in Plaquette 1. But due to the non-orthogonality of the plaquettes and the 
equidistance of the Cu-0, it can also form a triplet state with Plaquette 2 [19]. However, 
this is less likely because the ZRS hybridization is stronger and consequently more favoured 



14l |. Furthermore, the ZRS hybridization in plaquette 1 leads to spontaneous exchange of 



spins from the interaction part of Eq. (9). Therefore the effective hole seen by the Cu in 
Plaquette 2 has an opposite spin so that Eq. (9) become effective in forming new ZRS in 
this plaquette, resulting in a scenario in which the effective bound state forming the ZRS is 
propagated from one plaquette to another in the t-J model and thereby constitute a prob- 
able Cooper channel. This mode of formation and propagation of the ZRS is supported by 
the observation in [25| that the critical Tc of the superconducting cuprates depends on the 



number of Cu02 planes within a short distance of each other in the structure. 
Now the above description looks superficially rather too simplified to explain all the intrigu- 
ing properties of the high Tc in the suerconducting cuprates oweing to the present sophisti- 
cation of the literature which unknowingly may have become part of the complex nature of 
these materials. It is therefore necessary to emphasize that the mechanism proposed here do 
not rule out the effects of the other parameters and the influence of electron-phonon inter- 
actions from displacement of the atoms 26|]. For example, the spontaneous spin exchange 
between the O and Cu in one plaquette usually leads to virtual superexchange between the 
NN Cu atoms. Interestingly, such effects have been included as Cu-Cu exchange interactions 
and sometimes also as constituting additional hopping term (if:') in attempts to account for 
experimental data in several studies 11, [l^. Thus it is hoped that the formulation here 
gives the underlying physics of ZRS propagation in the Cu02 planes which can lead to more 
accurate investigation with methods that allow unbiased treatment of necessary parameters. 
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